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ABSTRACT 



We present deep V and I photometry of the isolated dwarf galaxy VV124=UGC4879, obtained from archival images taken with 
the Hubble Space Telescope - Advanced Camera for Surveys. In the color-magnitude diagrams of stars at distances larger than 
40" from the center of the galaxy, we clearly identify for the first time a well-populated old Horizontal Branch (HB). We show 
that the distribution of these stars is more extended than that of Red Clump stars. This implies that very old and metal poor 
populations becomes more and more dominant in the outskirts of VV124. We also identify a massive (M = 1.2 + 0.2 x 10 4 Mq) 
young (age= 250 ± 50 Myr) star cluster (CI), as well as another of younger age (C2, S 30 ± 10 Myr) with a mass similar to clas- 
sical open clusters (M < 3.3 + 0.5 x 10 3 Mq). Both clusters lie at projected distances smaller than 100 pc from the center of the galaxy. 
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1. Introduction 

The dwarf galaxy VV124=UGC4879 has been recently recog- 
nized to lie at the outer fringes of the Local Group (LG), near the 
turn-around radius (Ko pylovet al . 2008; Tikho nov et alJ uOlO). 
It is extremely isolated, being at =2 1 .3 Mpc from the nearest giant 
galaxy (the Milky Way) and at ^0.7 Mpc from the nearest dwarf 
(Leo AX 

In lBellazzini et al.1 (201 1, Pap-I hereafter) we presented new 
deep wide-field photometry of VV124, w hich allowed us to con - 
firm and refine the distance estimate by Kopylov et al. (2008), 
to study in some detail the stellar content of the galaxy and, 
in particular, its structure out to large radii, revealing the pres- 
ence of low surface brightness structures in the outskirts. We 
also presented the first detection of H 1 in VV124, obtained with 
the Westerbork Synthesis Radio Telescope (WSRT). 10 6 M 
of Hi were found, with a maximum column density of ^ 
50 x 10 19 cirT 2 . It is worth noting that the distribution of Hi 
is l ess extended t han th e stellar body of the galaxy. 

iJacobs et alJ d201ll Jll hereafter) presented Hubble Space 
Telescope (HST) - Advanced Camera for Surveys (ACS) ob- 
servations of VV124, resolving very well the stars down to 
F814W - 27 also in the most central region of the galaxy. 
In particular, they studied the stellar content in the innermost 
40" (see their Fig. 1), deriving the star formation history (SFH) 



Send offprint requests to: M. Bellazzini 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are asso- 
ciated with program GO-1 1584 [PI.: K. Chiboucas]. 



in that region by means of the synthetic co lor-magnitude di- 
amgram (CMP) techn ique (see, for example iRizzi et all 12002: 
Cignoni & Tosil l2010i and references therein). Jll concluded 
that most of the SF (^ 93 percent of the total, in mass) occurred 
at the earliest epochs, between 14 and 10 Gyr ago. After this, the 
SF ceased until 1 Gyr ago, when ^6 percent was formed between 
1 and 0.5 Gyr ago and nearly 1 percent in the last 0.5 Gyr. The 
derived SFH is similar to that observed in other nearby dwarf 
irregular/transition-type galaxies (see, for example IWeisz et al. I 
2011). The prevalence of old stars seems to increase at larger 
distances from the center (see also Pap-I). From the morphology 
of the blue-loop region, Jll inferred that the metallicity of the 
youngest stars is [M/H] - -0.7. 

In the present contribution, we use the same observational 
material of Jll to address three issues raised in Pap-I and not 
considered by Jl 1. These are, in particular: 

- The tentative detection of an extended Horizontal Branch 
(HB), including stars in the color range of RR Lyrae vari- 
ables and blue HB stars, in the CMD of the less-crowded 
outer regions of the galaxy (see Fig. 9 of Pap-I). While Jll 
report that a population of HB stars reaching V — I ~ 0.4 
is likely present in their CMD, its clean identification and 
quantification is hampered by the overlap with the young 
main sequence (MS) in the central region they consider. It 
would be interesting to search for the HB in the outer por- 
tion (R p > 40") of the ACS field, where young stars should 
be quite rare and the old component should be dominant. 

- A possible correlation was noted between the asymmetric 
distributions of the young stars and of the neutral Hydrogen. 
Since the young population is concentrated in the innermost 
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part of the galaxy, it can be traced in much deeper detail 
with the ACS data. This, in turn, prompts a more stringent 
comparison with the H i distribution. 
- We identified a candidate cluster located near the center of 
the galaxy and within the sheet of young stars (Cluster 1, CI 
hereafter). The integrated spectrum of the source, shown in 
Pap-I, shows strong H^ in absorption, suggesting a young 
age. With the high spatial resolution ACS image, we can 
hope to ascertain the real nature of CI and look for other 
conspicuous clusters, if any. 

There is no reason to repeat the excellent analysis of Jll, 
hence we focus strictly on the scientific cases listed above. 

In the following, we adopt the parameters of VV124 derived 
in Pap-I (see their Tab. 1), if not stated otherwise. In particu- 
lar, we adopt D - 1.3 Mpc (in good agreem ent with Jll) and 
E(B - V) = 0.015 (from ISchlegel et alj]l998h . At this distance 
1" corresponds to 6.3 pc, 1' to 378 pc. We also adopt the de- 
projected coordinate system X,Y of Pap-I, centered on the cen- 
ter of VV124 and rotated such to let the X axis coincide with the 
major axis of the galaxy. Finally, we use the elliptical radial dis- 
tance from the center of the galaxy r e . To give an example, stars 
fulfilling the condition r e < V are enclosed within an ellipse 
whose major axis coincides with the major axis of VV124 and 
has o = l' and e - 0.44 (see Fig.Q] below). The normal (pro- 
jected) radial distance from the center of the galaxy will be de- 
noted with R p . All theoretical models used for comparison with 
observed CM Ds are from the (canonical, solar- scaled) BASTI 
set dPietrinferni et al J 12004 iPercivalet all 12009). in the native 
ACS- Wide Field Channel (WFC) VEGAMAG systerrfl 

The plan of the paper is as follows: in Sect. [2] we describe the 
observational material and the data reduction process; in Sect. [3] 
we present the clear detection of the extended HB in the external 
regions of the ACS field and we discuss the radial distribution 
of the various stellar populations of VV124. Sect. [4] is devoted 
to a brief discussion of the correlation between the spatial distri- 
bution of the H i and of the young stars, and in Sect. [5]we study 
the candidate luminous star clusters identified here and in Pap- 
I. Finally we briefly summarize and discuss our main results in 
Sect.|6] 




Fig. 1. Upper panel: stars detected in the ACS images (dark dots) 
are superposed to the stars from the LBT sample of Pap-I (grey 
dots) in a projected coordinate system centered on the center 
of VV124 and rotated to have the X axis coinciding with the 
major axis of the galaxy (see Pap-I). Only stars brighter than 
/ = 26.5 from the ACS sample have been plotted, for clarity. 
Lower panel: zoomed view of the central part of the ACS sample 
in the same coordinate system. Two concentric ellipses of semi- 
major axis a — 1 .0' and a - 1 .5', having the same orientation and 
ellipticity of the galaxy isophotes (see Pap-I) are overplotted (in 
white and in black, respectively). The white circle has a radius 
R p - 40", enclosing the region studied by Jll. Blue stars lying 
on the young Main Sequence (-0.4 < V - I < 0.0 and / < 25.5, 
see Fig. [3} are plotted as red filled circles, the brightest ones 
(having / < 23.0) are encircled in orange. 



2. Observations and data reduction 

The observational material consists of two pairs of F606W and 
F814W ACS-WFC0 images obtained within the program GO- 
1 1584, which we retrieved from the HST archive. The total ex- 
posure time of the stacked drizzled images is 1 140 s and 1 150 s 
for the F606W and F814W images, respectively. The WFC is a 
mosaic of two 4096 X 2048 px 2 CCDs, with a pixel scale of 0.05 
arcsec px" 1 . The center of VV124 is placed at the center of one 
of the mosaic chips (chip 1). 

Positions and photometry of individual stars were obtained 
with the ACS module of the Poi nt Spread Funct ion (PSF) - fit- 
ting package D OLPHOT v.l .lB (Dolphin 2000), as described, 
for example, in iGalleti et al.l (2006). The package identifies the 
sources on a reference image (in this case the drizzled F814W 
image, as in Macke v et alJl20 06) and performs the photometry 
on individual frames, also taking into account all the information 
about image cosmetics and cosmic-ray hits which is attached to 
the observational material. We adopted a threshold of 3cr above 



1 http://193.204. 1.62/index. html 

2 http://www.stsci.edu/hst/acs/ 

3 http : //purcell . as . arizona . edu/dolphot/ 



the background for the source detection. DOLPHOT automat 
ically provides magnitudes both in the WFC VEGAMAG sys 
tem and in the Johnson-Kron-C ousins system (Dolphin 200< 
As done in IGalleti et al.l d2006l) and iPerina et alj ( 120091 I20T 
we selected the sources according to the quality parameters pro- 
vided by DOLPHOT, to clean the sample from spurious or badly 
measured entries. In particular, we retained in the final cata- 
logue only stars having quality flag-1 (i.e. best-measured stars_ 
chi < 2.0, |sharp| < 0.5, crowding parameter crowd < 
and photometric error < 0.3 mag in both filters (see Dolphin 
2000, for details on the meaning of the various parameters). In 
the following sections we will always adopt this selected cata- 
logue that contains 44894 sources, in total 0. 

We transformed the star positions from this system into 
J2000 RA and Dec by fitting a fourth degree polynomial to 5038 
stars in common between the ACS catalogue and the LBT cata- 



4 We will relax this selection criterion to crowd < 1.0 in Sect. 5, to 
study in more detail the star clusters, where stars are significantly more 
crowded than average. 

5 This catalogue is available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/.../ .. 
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Fig. 3. Upper panels: CMD of VV124 from the ACS data within (left panel) and outside (right panel) a circle of radius R p - 40". 
The average photometric errors are reported as error bars in the left panel. Lower panels: CMDs in other regions of the ACS field 
delimited by the ellipses plotted in Fig. Q] The right panel regards a radial region similar to that considered in Pap-I to tentatively 
detect the extended HB population (see their Fig. 9). In all the panels stars are plotted as black points in regions of the CMD with 
few stars and as grey squares otherwise, with the scale of grey proportional to the local density of stars in the CMD (lighter tones of 
grey correspond to higher density). 



logue of Pap-I, using the CataPack suite of codeo The r.m.s. of 
the residual was 0.07" in both RA and Dec. The positions of the 
stars in the ACS field with respect to the LBT sample are shown 
in the upper panel of Fig. Q] In the lower panel of the same figure 
we zoom on the part of the ACS field containing the center of the 



6 Developed by Paolo Montegriffo at INAF 
Astronomico di Bologna. 
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galaxy. It is particularly interesting to note that a large number 
of stars are detected beyond the R p < 40" limit of the sample 
used by Jl 1 for their determination of the SFH. It is beyond that 
radius, and especially beyond r c = 1.5' where we expect to de- 
tect the extended HB of VV124, an obvious tracer of very old 
(ageS; 10 Gyr) and metal-poor populations. 
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Fig. 2. Results of the artificial stars experiments. Upper panels: 
Difference between the output and input magnitudes of the ar- 
tificial stars as a function of F814W magnitude. Lower panels: 
Completeness as a function of F606W magnitude in the inner 
(r < 1.0'; left panel) and outer (r > 1.0'; right panel) regions 
of the galaxy. Red continuous curves shows the completeness 
for stars in the color range of the RGB and the RC, blue dashed 
curves refer to stars in the color range of the old HB and the 
YMS (see Fig. El below). The blue curves start at F6Q6W = 25.0 
because in that color range we produced only artificial stars 
fainter than this limit, focussing to the study of the complete- 
ness in the blue HB region. 

2. 1 . Artificial stars experiments 

To study the completeness of the sample and the photometric 
uncertainties, we performed a set of arti ficial star experim ents, 
following the same procedure adopted bv lPerina e t al. (200§l see 
that paper for a detailed description and references). We pro- 
duced =* 200000 artificial stars having 21.0 < F6Q6W < 28.0, 
with colors in the range covered by the bulk of observed stars, 
i.e. -0.4 < F606W - F8UW < 1.3. 

The results of the experiments are presented in Fig. El The 
upper panels show that the photometric accuracy is better than 
0.1 mag for F6Q6W S 25.0; a quantitative summary of the re- 
sults displayed in these panels is presented in TableQ] The lower 
panels show the completeness factor (C/) as a function of mag- 
nitude for two radial ranges and for two color ranges. At mag- 
nitudes fainter than F606W - 25.0, the completeness is slightly 
worse for blue stars than for red stars. However, even blue stars 
in the inner radial bin have Cf ^50% for F6Q6W ;£ 26.5. In 
the outer radial bin, the same level of completeness is reached at 
F606W =* 27.2. 



3. The extended HB of VV124 

In the upper left panel of Fig. El we show the ACS CMD of 
VV124 in the same radial range studied by Jl 1 (R p < 40"). The 
two CMDs are very similar, as expected, since they are obtained 
from the same dataset and with the same photometry package. 



Table 1. Photometric errors from artificial stars experiments 



F814W cr(FSUW - FSUWt) <r{col - colt)" 



20.0 6 


0.009 


0.006 


21.0 


0.010 


0.009 


22.0 


0.013 


0.014 


23.0 


0.021 


0.024 


24.0 


0.040(0.045) c 


0.044(0.045) 


25.0 


0.077(0.085) 


0.084(0.087) 


26.0 


0.153(0.164) 


0.175(0.180) 


27.0 


0.236(0.242) 


0.281(0.286) 



Notes. ° col = F606W o - F814W ; coh = F606W, - FUAW t . 
b The reported errors are standard deviations computed for recovered 
artificial stars having F814W„ within ±0.5 mag from the value reported 
in column 1 . 

c Numbers within parentheses refer to the sample with CRO < 1.0, 
which is adopted in the analysis of the clusters. All the other values 
refer to the standard CRO < 0.3 sample. The alternative value is re- 
ported only in the cases when it differs by more than 0.001 mag from 
the standard value. 



The most obvious features are the broad RGB, with the RGB 
tip located at F814-W ~ 21.5, and the Red Clump of intermedi- 
ate/low mass core He burning stars around F814W ~ 25.3. A 
handful of Asymptotic Giant Branch stars is present on the ex- 
tension of the RGB above the tip. The obvious blue plume of 
young MS stars runs over all the considered magnitude range 
around F606W - F814-W ~ -0.2. The sparse plume protrud- 
ing from the RGB at F814W ~ 23.5 toward the blue has been 
identified by Jl 1 as composed of Blue Loop stars of metallicity 
[Fe/H]~ -0.7; these are the evolved counterpart of the young 
MS. 

Jl 1 noted that the synthetic CMD best reproducing the ob- 
servations included HB stars reaching colors as blue as V - I ~ 
F606W - F814W =* 0.4 that were not identified in the ob- 
served CMD. We note, on the other hand, that it seems possi- 
ble to discern a blue extension of the HB in our CMD around 
F&UW ~ 25.8 and 0.0 < F606W - F814W < 0.6. This impres- 
sion is confirmed by the fact that the CMD of stars beyond the 
region studied by Jl 1 display an obvious extended HB (upper 
right panel of Fig. E), with an higher-density portion between 
F606W - F8 UW =* 0.6 and F6Q6W - F8 UW - 0.3, and a less 
populated part reaching F606W - F8UW =* -0.4. 

The extended HB becomes more evident when more external 
regions are considered, as shown in the lower panels of Fig. El 
This should be partially due to the decreased crowding, which al- 
lows to get more accurate photometry, thus decreasing the width 
of the other CMD sequences that may partially overlap with the 
HB. However, the most important effect is the decreasing impact 
of the young MS plume, which contributes to hide the blue HB 
near the center of VV124. We will show below that, in addition 
to these observational factors, there is also a genuine population 
gradient making the HB stars more abundant in the outer regions 
of the galaxy. The lower right panel of the figure shows the CMD 
in the same radial region where we obtained the tentative detec- 
tion of the extended HB from the LBT data, thus providing a 
direct confirmation of that finding. 

In the upper right panel of Fig. |4] we superposed a synthetic 
HB population, produc ed from the BASTI theoretical models 
(Pietrinfer ni et al.ll2004l) with the dedicated web tooQ, to the ob- 
served CMD of the outermost region of VV124 sampled by the 



7 http://193.2Q4. 1.62/index. html 
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Fig. 4. Upper left panel: The adopted selection boxes are plot- 
ted on the CMD for r e > 0.5' (the innermost region has been 
excluded to avoid confusion in the panel). Upper right panel: 
A synthetic HB population computed from BASTI models, re- 
ported to the same distance and reddening of VV124, is super- 
posed to the observed CMD for r e > 1.5'. The radial selec- 
tion is adopted to obtain the cleanest view of the HB. Red HB, 
RR Lyrae variables, and Blue HB are plotted in red, green and 
blue, respectively. Lower panel: Radial distribution of the vari- 
ous populations selected with the CMD boxes, as indicated by 
the labels. 



ACS data. The synthetic population (1000 stars) have Z = 0.001 
(appropriate for the old population of the galaxy, according to 
Kopylov et al. 2008, Pap-I and J 1 1), a mean mass of 0.494 Mq 
and a dispersion cr = 0.2 Mq, and it has been shifted to the 
distance and reddening of VV124. This is not intended to re- 
produce the observed HB but only to provide a sanity check for 
the interpretation of the sequence and to have a hint on the clas- 
sification of stars within the HB (color coded according to the 
classification provided by the model). In particular, it is evident 
that the galaxy likely host a significant population of RR Lyrae 
variables, since there are many HB stars in the color range typ- 
ical of these variables. Clearly their definitive classification as 
genuine RR Lyrae would need time series photometry, which is 
not available to us. The detection of such a rich population of 
candidate RR Lyrae and blue HB stars is fully compatible with 
the SFH derived by Jl 1, who concluded that <: 90% of the stars 
were formed between 14 and 10 Gyr ago. 

For F606W-F814W^ 0.0, the bluest portion of the ob- 
served HB seems to depart from the behavior predicted by 
the models, showing only a modest decline of the mean mag- 
nitude with color, while the synthetic population plunges to- 
wards much fainter luminosity. This may be due to a combi- 
nation of effects, including non-perfect temperature-color trans- 
formations in this extreme color regime and photometric er- 
rors/completeness/blending/selection effects due to the proxim- 
ity of the limiting magnitude level, which becomes brighter for 
bluer colors. In Fig. [5] we simulate the observational effects on 
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Fig. 5. Left panel: CMD of the synthetic stars from the BASTI 
model shown in Fig.|4]having 25.5 < F6Q6W < 28.0. The mean- 
ing of the symbols is the same as in Fig. [4] The continuous line is 
the corresponding Zero Age HB model, plotted in both panels as 
a reference. Right panel: simulation of the observation process 
on the synthetic stars, obtained from the artificial stars experi- 
ments (sample with r e > 1.5'). Note that the output HB mor- 
phology is very similar to the observed one. 



a subsample of the synthetic HB stars. The experiment demon- 
strates that the apparent mismatch can be fully accounted for by 
the effects of photometric errors and blending. Hence, most of 
the stars in the horizontal strip around F6Q6W ~ 26.5 and bluer 
than F606W- F814W ~ 0.0 are genuine BHB stars (at least for 
r e > 1.5', where the contamination by YMS stars is negligible), 
in spite of the apparent mismatch with the raw modelo 

The upper left panel of Fig. [4] shows the boxes we adopted 
to select stars on the Young Main Sequence (YMS), in the Red 
Clump (RC) and in the extended HB, to compare their spatial 
distributions. For the HB we adopted a box covering a rela- 
tively narrow color range, approximately matching the one of 
RR Lyrae variables. This choice was forced by the need to avoid 
contamination from RC stars, on the red side, and, above all, 
from faint YMS stars on the blue side. This would have been par- 
ticularly problematic since (a) YMS are obviously much more 
centrally concentrated than older stars, and (b) the YMS blue 
plume becomes wider and reaches redder colors at faint mag- 
nitudes (V ^ 25), clearly overlapping with the blue HB in the 
central region. We verified that the actual radial distribution of 
HB stars bluer than the adopted selection box is fully compatible 
with a population having the same radial distribution as the HB 
sample plus a strong contamination of YMS in the innermost 
~ 0.5'. Hence the selected HB sample can be considered as a 
fair tracer of the old HB population as a whole, including BHB. 
The selected RC and HB samples should have similar complete- 
ness properties, since they have similar magnitudes and colors. 
To verify this hypothesis we compared the behavior of the com- 



8 To simulate the observations we proceeded as follows. For each 
model star plotted in the left panel of Fig. \5\ we isolated the subset of 
artificial stars having input F606W magnitude within ±0.2 and input 
color within ±0.05 of the considered synthetic star. Then we extracted 
at random one of the selected artificial stars. If the extracted artificial 
star is not recovered in the output we remove the corresponding syn- 
thetic star, thus mimicking the effect of incompleteness. If it is recov- 
ered, we compute its output-input difference in magnitude and color 
and we sum these quantities to the magnitude and color of the syn- 
thetic star, thus mimicking the effects of photometric errors and blend- 
ing (see Bellazzini et al. 2002b, for a similar application). We repeated 
the experiment several times and we verified that the HB morphology 
in output does not depend on the specific random realization. For an 
example of detailed modeling of an HB population in a similar context 
see lMonelli et ali ( 120 id . and references therein). 
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pleteness as a function of r e for the two samples, using our set 
of artificial stars. It turned out that the two completeness curves 
are undistinguishable: hence any observed difference in the ra- 
dial distribution of RC and HB stars cannot be accounted for by 
completeness effects. 

The cumulative radial distributions of the three selected pop- 
ulations are compared in the lower panel of Fig. [4] In addition 
to the striking (and already noted) difference between the YMS 
and the older RC and HB populations, it is clearly revealed that 
the HB stars have a more extended distribution with respect to 
RC stars. A Kolmogorov-Smirnov test gives that the probability 
that the two samples are extracted from the same parent popu- 
lation is P < xlCT 10 , thus firmly establishing the significance 
of the difference. It must be concluded that VV124 displays 
the classical population gradient observed in all dwarf galaxies, 
with the older/more-metal-poor populations becoming increas- 
ingly more abundant at larger radii and youn ger/more-metal-rich 
being more concent rated toward the center (Harbeck et al. 2001; 
IBellazzini etai]|2002al l200l iTolstov. Hill & Tosi et al.ll2009h . 

3.1. The color distribution of RGB stars 

In Pap-I, the comparison of the observed RGB with empirical 
templates from Galactic globular clusters revealed the possible 
presence of a small metallicity gradient in the same sense. The 
analogous comparison, performed with the much more accu- 
rate ACS photometry is presented in Fig. [6] to check the find- 
ings of Pap-I and to verify whether the distribution and mean 
color of RGB stars is compatible with the gradient found with 
the HB stars. It must be stressed that the following analysis 
is not intended to provide accurate metallicity estimates, since 
it is well known that photometric determinations of the metal- 
licity can be affected by serious degeneracies when the under- 
lying age distribution is not well constrained (see ICole et alTI 
2005; Lianou, Grebel & Koch 2011, and references therein). 
On the other hand we look for differences in the inner and 
outer samples of RGB stars that may or may not be in agree- 
ment with the observed HB gradient. In the analysis we will 
adopt the global metallicity parameter [M/H] introduced by 
Salaris, Chieffi & Straniero (1993) to allow meaningful com- 
parisons between the RGB colors of stell ar populations hav- 
ing different g-element abundances (see also Ferraro et all 1999: 



IBellazzini et~al1l2002atlSalaris. Cassisi & Weissll2002l) . 

Fig. [6] shows that the fraction of stars spilling to the red of 
the [M/H] =-1.21 template is indeedlargerin the 0.5' < r € < 1.0' 
range than at larger distances, even if the amplitude of the color 
spread is smaller than that apparent in Pap-I. Considering only 
stars in the range 2 1 .7 < F814W < 23.0 (i.e. the least affected by 
photometric errors and crowding effects) and enclosed between 
the bluest and the reddest template, we computed the fraction 
of stars lying to the red of the [M/H] =-1.21 ridge line. This is 
0.29 + 0.03 in the range 0.5' < r e < 1.0' and 0.20 + 0.02 for 
r e > 1.0' (the same result is obtained also if stars with r e < 0.5' 
are included in the inner bin, and if stars down to F% 14W = 23.5 
are considered). This is quite similar to the results obtained in 
Pap-I and in broad agreement with the gradient shown in Fig. [4] 

We obtained metallicity estimates for individual stars with 
F814W< 23.0, by linear interpolation within the templates grid, 
as in Pap-I. The resulting metallicity distributions are quite 
broad, ranging from the blue limit of the grid, [M/H] =-1.95, to 

[M/H] 1.0. The mean metallicity in the range 0.5' < r e < 1.0' 

is ([M/H]) = -1.4, while for r e > 1.0' it is ([M/H]) = -1.5. 
The associated standard deviation, deconvolved from the obser- 
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Fig. 6. Comparison of the CMDs of VV124 in different (ellipti- 
cal) radial annuli with the RGB ridge lines for the Galactic glob- 
ular clusters NGC 63 41, NGC 675 2 and NGC 104 (from left to 
right), from the set by B edin et alj (2005). that provides the tem- 
plates in the same ACS-WFC VEGAMAG system as our obser- 
vations. The reddening, distance moduli and metallicities of the 
clusters are taken from Ferraro et al. (1999), as in Pap-I. In the 
right panel the ridge lines are labeled with the [M/H] va lue of 
the corresponding c luster in the Carretta & Grattorj (Il997l) scale 
(col. 4 of Tab. 2 of lFerraro et alj|1999l) . The radial range of the 
left panel coincides with the innermost radial range considered 
in Pap-I. 



vational effects, is o-[m/h] — 0.25 dex in both radial bins, also in 
reasonable agreement with the results of Pap-I. 

4. The spatial distributions of H 1 and young stars 

In Figure [7] we compare the location of those stars identified 
to be on the YMS with the Hi distribution. This comparison 
clearly shows that, when considered on large scales, there is a 
very good spatial correspondence between the distribution of the 
YMS stars and that of the H 1 , both show the same asymmetry 
with respect to the overall light distribution of VV124. Such an 
overall correspondence is not entirely unexpected, because, af- 
ter all, stars form from gas. On smaller scales (i.e. less than 100 
pc) the spatial correspondence between YMS stars and the Hi is 
less good, with the peak of the H 1 somewhat displaced from the 
YMS stars. This situation is not untypical for small galaxies: a 
good correlation between star formation and H 1 on large scales 
and the breakdown of it on smaller scales is observed in many 
small galaxies (Begum et al. 2006; Rovchowdhury et al. 2009). 
The close overall correspondence between the YMS and the 
H 1 also suggests that the observed kinematics of the H 1 is also 
connected to the formation of the stars we now observe to be 
on the YMS. As discussed in Pap-I, if the Hi tail is interpreted 
as an outflow, the energy output of roughly a dozen OB stars 
would be sufficient to drive it. However, an interesting alterna- 
tive is that the tail corresponds to infall of gas and that VV124 
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is re-accreting part of the gas which was expelled by the for- 
mation of the stars we now see as an old population. That such 
re-ac cretion may occur is suggested by some recent simulations 
(e.g. lSawala et al.|[201 lb . and would also be very well consistent 
with the fact that VV124 has been dormant for a long period of 
time since the early phase of star formation and only very re- 
cently started to form young stars again. The metallicity of the 

young stars in VV124 ([Fe/H] 0.7, as estimated by Jll from 

the observed morphology of the Blue Loop) also suggest that, 
if they formed from the HI gas cloud, this cloud must contain 
enriched material and hence it consists of re-accreted material. 



5. Luminous star clusters in VV124 

The zoomed view of the galaxy center shown in Fig. [8] shows 
that the candidate star cluster identified in Pap-I (Cluster 1, CI 
hereafter) is indeed largely resolved into stars in the ACS im- 
ages, hence it is likely a genuine star cluster (see below). On 
the other hand, the two bright sources to the west of CI are a 
background galaxy and a foreground star, respectively. There are 
several small groups of young stars in this regions that may be 
possibly classified as open clusters or associations. However, the 
only other group of size/luminosity comparable to CI is the one 
indicated as Cluster 2 (C2) in the figure. C2 appe ars to coincide 
with th e position of the H 11 region identified by Ko pylov et al.l 
(2008, by means of a long slit optical spectrum). With the avail- 
able data it is not possible to establish whether the H 11 region is 
associated to the cluster or it is just a chance superposition. 

In Sect. 5.1 and 5.2, below, we will study in deeper detail 
these two small stellar systems. The insight into the physical 
properties of the clusters is greatly hampered by stellar crowd- 
ing which is maximized with respect to the surrounding field. 
Indeed, to reach F814W ~ 26.5 in the very small areas en- 
closing the clusters, we had to relax the selection on the crowd- 
ing parameter, accepting stars up to crowd < 1.0. In addition, 
the contamination by field stars provides additional problems 
to the interpretation of the CMDs. For these reasons, the best 
we can hope to obtain from the analysis of the clusters CMD 
are age estimates with uncertainties of order 20%-50%, also 
due to the uncertainty in the associated metallicity and the sys- 
tematics associated with uncertainties in the theoretical models 



(Gal lart. Zoccali & Apariciol l2005). However such large age un- 
certainties translates into uncertainties of 10%-20% in the mass 
estimates obtained by converting total luminosity into mass (see 
iPerina et al . 2009, 2010, and references therein). 

In the following analysis we will try to estimate ages by 
fitting the cluster CMD with BASTI isochrones, following 
IPerina et all Q2P09). As a reference value for the cluster metal- 
licity we assume Z = 0.004, i.e. the mean value estimated by 
Jll for the young population of VV124. Since the two clus- 
ters are clearly associated with this young component, this is 
by far the best guess on the metallicity of the clusters. While we 
will present only solutions with Z = 0.004, we have attempted 
to fit the observed CMDs also with isochrones with metallicity 
ranging from Z = 0.001 to Z = 0.02 and we found that the 
Z = 0.004 models always provide the best overall fit. Estimates 
obtained with the Girardi et al. (2000) isochrones, for compari- 
son, gives age values larger by ~30% that those obtained with 
BASTI isochrones (hence slightly larger mass estimates), keep- 
ing fixed any other parameter. Finally, the good match in color 
between the MS of the models and of the clusters shows that the 
reddening value adopted for the whole galaxy is appropriate also 
for the clusters, indicating they are not affected by additional ex- 
tinction. 



5.1. Cluster 1 

In Fig. we compare the CMD in a small area (inner 
radius=0.6", outerradius=2.5") centered on CI with those from 
identical areas centered on positions symmetric with respect to 
the center of V V 1 24 and which are used as control fields (CF, see 
the upper right panel of the figure). It is quite clear that, even if 
the CI region has a larger area devoid of measured stars because 
of the higher crowding associated to the cluster center, the CI 
CMD show an excess of MS stars fainter than F814W - 24.5 
and, in particular, an obvious excess of BL stars, with respect 
to all the considered CFs. This is clearly observed even if the 
analysis is limited to stars with crowd < 0.3. In particular in 
the CI CMD there are 18 stars lying in the superposed BL se- 
lection box, while 3, 5, and 4 are found in the CFg, CF^, and 
CF\y CMDs, respectively. The associated excess of BL stars in 
CI with respect to the three CFs is at the 3.3<x, 2.7cr, and 3.0cr 
level. If we limit to the sample of stars with crowd < 0.3, we 
find 1 1 stars in the BL box of the CI CMD, while only 3 stars 
are found in the CMD of each CF. 

In Table|2]the ratio of stars bluer than F606W-F814W = 0.2 
(N B ) to the stars redder than F606W - F8 14W = 0.4 (N R ) in the 
range 24.0 < F814W < 26.0 is shown to be significantly larger 
(at a > 3<r level) in CI than in any considered CF. Given the 
regular, compact and roundish appearance of CI and the obvious 
excess of young stars found in the CI area, we conclude that 
CI is a genuine star cluster. It is worth noting that the BL and 
the MS fainter than F814W ~ 24.0 are the only features of the 
CI CMD that can be unequivocally attributed to the cluster, as 
they are significantly over-abundant in the CI field than in the 
CFs. On the other hand, one can be tempted to interpret the five 
stars of the CI CMD lying in the region F6Q6W - F814W < 
0.2 and F814W < 24.0 as very young MS stars associated to 
the cluster. However this interpretation has no statistical basis, 
since in the same region there are 3, 1, and 3 stars in CF£,CFn, 
and CFw, respectively. Hence this handful of bright and blue 
stars in the CI CMD is fully compatible with the expectations 
for the field population. Moreover, we did not find an isochrone 
that simultaneously matches these very blue stars and the BL 
population, that is unequivocally associated with CI. 
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Fig. 8. The two candidate clusters (CI and C2) are indicated on this image of the center of VV124 (from the drizzled F606W image). 
The circle around CI has a radius of 2.5", the one around C2 a radius of 2.0". The ellipse has the same ellipticity and orientation of 
the VV124 isophotes and semi-major axis a = 25". The center of the galaxy is indicated by a cross. 

Table 2. %r ratios for the candidate star clusters and their control fields 



Cluster Name 


Cluster 


CF £ 


CF W 


CF W 


Min. Excess c 


Max. Excess d 


CI" 

C2 6 


0.84±0.16 
0.70±0.16 


0.11 ±0.04 
0.22±0.06 


0.13±0.04 
0.13±0.04 


0.28± 0.06 
0.29± 0.08 


3.3tr 
2.3tr 


4.4o- 
3.4cr 



Notes. " N B is the number of stars having F606W - F814W < 0.2 and 24.0 < FS14W < 26.0, while N R is the number of stars having F606W - 

F814W > 0.4 in the same magnitude range. 

b The definitions of N B and N R are the same as for CI, except for the magnitude range that is 22.0 < F814W < 26.0, in this case. 

c Difference between the jf- in the cluster and in the CF showing the lowest value of this parameter, in units of cr. 

d Difference between the jf- in the cluster and in the CF showing the highest value of this parameter, in units of cr. 



In the upper left panel of Fig. [9] we show also the compari- 
son of the observed CMD of CI with the isochrone correspond- 
ing to our preferred solution. In the age range covered by the 
young population of VV124, the luminosity of the Blue Loop is 
a very sensitive age indicator (Dohm-Palmer et al.ll997l) : the BL 
of C 1 is well matched by the adopted isochrone, having an age of 
250 Myr. Such a young age is in agreement with the integrated 
spectra of CI we presented in Pap-I. Following Perina et al. 
(2009), we associate an uncertainty of +50 Myr to our age es- 
timate. All age estimates we obtained by adopting isochrones of 
different metallicity (in the range Z=0.001 -0.02) lie within this 
error-bar. 



responding to 15.7 pcQ . These, as well as other relevant cluster 
parameters, are listed in Tab. [3] While the cluster may be slightly 
more extended than the adopted radius, this choice should be 
a reasonable compromise between including most of the clus- 
ter light and avoiding strong contamination from the dense field 
population. The sky level for aperture photometry was estimated 
in a nearby region, to subtract also the contribution from field 
stars of VV124. From aperture photometry on concentric annuli, 
and assuming r - 3.0" as the limiting radius of the cluster, we 
also obtained a rough estimate of the half-light radius, r/, =* 1.5", 
corresponding to =* 9.5 pc, has been consistently obtained from 
both the F606W and the F814W images. This is within the range 



We obtained estimates of the integrated V and I magnitudes 9 Integrated magnitudes have be en transformed from F606W,F814W 
from aperture photometry within a circular radius r = 3.0", cor- to V,I using the transformations by Siria nni et"aD (2005) 
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Table 3. Observed and derived parameters of cluster CI 



Parameter 


value 


notes 


a 


09:16:03.14 


J2000* 


So 


+52:50:27.0 


J2000* 


age 


250 ± 50 Myr 


a 


[M/H] 


-0.7 


a 


R P 


4.3" 


b 


''/, 


~ 1.5" 


c 


v,„, 


18.4 ±0.1 


d 


hoi 


18.1 ±0.1 


d 


My 


-7.2 ±0.16 




u 


6.5 ± 1.0x10*1© 


total V luminosity 


*'*s1ars 


1.2 ± 0.2 x 10 4 M o 


stellar mass'' 



Notes. * In the astrometric system embedded into the F606W drizzled 
fits image. " From the isochrone fit of Fig. [9] (See Sect. 5). * Projected 
distance from the center of VV124. c From aperture photometry on con- 
centric annuli, assuming a limiting radius Ti = 3.0". d From aperture 
photometry with radius r = 3.0". e Assuming (M/L) v = 0.188 from 
the BASTI database, for a model having Age=250 Myr, Z=0.004 and a 
Kroupa IMF. 



12 12 12 

F606W-F814W F606W-FB14W F606W-FB14W 



Fig. 9. Comparison of the CMDs of the circular annuli around 
Cluster 1 and three nearby Control Fields indicated in the up- 
per right panel. In this map the x sign indicates the center 
of the galaxy. The plotted points correspond to sources hav- 
ing crowd < 1.0. In the CMDs, sources with crowd < 0.3 
are marked by an empty circle concentric to the filled one. An 
isochrone providing a reasonable fit to the observed CMD is also 
over-plotted to the CI CMD (thick red line). A selection box for 
the BL population is superimposed to each CMD. 



of rh values measured bv lBarmbv et alJ J2009I) in M31 clusters 
of similar age and mass (see below). 

The total V luminosity is L v = 6.5 + 1.0 x 10 4 Lq. 
Assuming ^ = 0.188, from a Simple Stellar Population (SSP 
iRenzini & Fusi Peccil fT988b model having Z = 0.004, age = 
250 Myr, we obtain Af» = 1.2 + 0.2 x 10 4 MgH This implies 
that CI is significantly brighter than Galactic open clusters of 
similar age and it is more massive than Galactic open clusters of 
any age (Bellazzini et al. 2008b; Perina et al. 2009), except for 
v ery few (and much younger) young massive clusters discussed 
in lFigerl d2008) and Messineo et al.l d2009t) . On the other hand, it 
appears very similar in age, mass and size to the young massive 
clusters found in relatively large numbers in many star forming 
galaxies of any mass (see Larsen & Richtler 1999; Perina et al. 
2010; iBarmbv et al.ll2009l: IVanzi & Sauvage 2006; Larsen et al. 
201 ll for references and discussion). In t he clu ster c lassification 
schem e adopted by lLarsen & Richtlerl d2000t) and iBillett et a] 
(2002) CI clearly belongs to the class of regular star clusters 



:ly 

rA 



(M v < -8.5, see Annibali et al. 2009). 



of young stars in C2 with respect to all the CFs, albeit at a lower 
level of significance than in the case of CI (Tab|2]i. We conclude 
that also C2 is likely a genuine cluster. 

It is very hard to age-date this cluster. If the H n region is 
indeed physically associated with C2, the age is constrained to 
be^ 10 Myr, and the three stars brighter than F814W = 21.0 are 
not cluster BL stars. This is clearly possible since, even if two of 
them are redder than F606W - F814W = 0.0 and do not have 
direct counterpart in the CMD of the CFs, this may be merely 
due to fluctuations inherent with such small number of stars. 

On the other hand, if these stars are interpreted as genuine 
cluster BL stars, the age should be larger than the above value 
and the coincidence with the Hn region must be considered 
as due to chance superposition. Unfortunately, the accuracy of 
the positioning of the H n region is not sufficient to check the 
degree of spatial co incidence with C2 o n scales smaller than 
1" -2" (see Fig. 2 in lKopvlov et al . 2008). Considering the three 
brightest stars as cluster members, the best fit is obtained with 
an isochrone of age 30+10 Myr. Given the above discussion, we 
take this value as an upper limit to the age of C2. 

The total V luminosity is L v - 5.4 + 0.8 x 10 4 Lq, the 
half-light radius is r;, =* 0.6", corresponding to ~ 4 pc (all 
the estimates as been obtained in the same way as for CI, 
see above). Assuming j^- - 0.061, from a SSP model hav- 
ing Z=0.004, age=30 Myr from the BASTI set , we obtain 
M, < 3.3 + 0.7 x 10 3 Mq. This is in the high-mass side of the 
Galactic open clusters distribution, but not out of the range cov- 
ered by those stellar syst ems. Also C2 should be included amon g 
regular star clusters dLarsen & Richtlerl2000llBillett et al2 002). 



5.2. Cluster 2 

C2 appears more irregular in shape with respect to CI, being 
clearly dominated by a few bright stars. This suggests a younger 
age than that of CI. We performed the same kind of analysis as 
done for the latter cluster. The jf- ratio indicates an overdensity 



10 All the BASTI SPSS models adopt a lKroupd d200lD initial mass 
function (IMF). 



6. Summary and Conclusions 

We have presented a new reduction and analysis of archival HST- 
ACS images of the isolated dwarf galaxy VV124. A detailed 
study of the SFH in the innermost regions of the galaxy, from the 
same data, has been already presented by another group (Jll). 
Here we have focused our attention on a few specific issues that 
have been raised in Pap-I and were not addressed in the Jll pa- 
per. 
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Table 4. Observed and derived parameters of cluster C2 



Parameter 


value 


notes 


a 


09:16:03.69 


J2000* 


So 


+52:50:20.4 


J2000* 


age 


< 30 ± 10 Myr 


a 


[M/H] 


-0.7 


a 


R ? 


11.9" 


b 


n, 


~ 0.6" 


c 


v,„, 


18.6 + 0.1 


d 


'tot 


18.6 + 0.1 


d 


My 


-7.0 + 0.16 




Ly 


5.4±0.8xl0 4 L© 


total V luminosity 


"1 stars 


< 3.3 ± 0.5 x 10 3 M Q 


stellar mass 1 " 



Notes. * In the astrometric system embedded into the F606W drizzled 
fits image. " From the isochrone fit of Fig. |9](See Sect. 5). * Projected 
distance from the center of VV124. c From aperture photometry on con- 
centric annuli, assuming a limiting radius r\ = 2.0". d From aperture 
photometry with radius r = 2.0". e Assuming (M/L) v = 0.061 from 
the BASTI database, for a model having Age=30 Myr, Z=0.004 and a 
Kroupa IMF. 
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Fig. 10. The same as Fig. [9] for Cluster 2. In this case the CFjy 
field is not symmetric with respect to the other CFs, to avoid 
inclusion of a significant portion of CI, whose position is also 
indicated in the map. 



- We provide a clear detection of an extended Horizontal 
Branch population. It has been shown that HB stars have a 
spatial distribution which is less concentrated than the RC 
population. This strongly supports the idea that very old and 
metal poor stars become more and more dominant at larger 
distances from the center of VV124. 

- The near coincidence of the asymmetric sheet of young stars 
near the center of VV124 and the density peak of the (also 
asymmetric) associated Hi cloud strongly support the idea 
that the complex kinematics and velocity gradients observed 
in the gas in Pap-I are most likely associated with the dy- 
namics of the star formation episode. It is interesting to 



note that some of t he isolated dwarfs in the simulations by 
ISawala et al.1 (120111) lose all their H i at very early times and 
are able to re-accrete some of the expelled gas at later times. 
This scenario would be fully consistent with the observed 
properties of the Hi in VV124 and would provide a natural 
explanation for the long period of quiescence between the 
early phase of star formation and the very recent sprout of 
activity that originated the sparse population of young MS 
stars currently observed (see Jl 1). 
- We confirm that the candidate cluster identified i n Pap-I (CI) 
is a g enuine young massive star cluster (see iPerina et all 
2010, and references therein). It has an age of ~ 250 Myr and 
a total mass of =* 1 .2 x 10 4 Mq, The cases of such low-mass 
early-type dwarf galaxies with massive clusters are quite rare 
in the Local Group (M98). We identified another likely clus- 
ter (C2), significantly younger (age< 30 Myr) and less mas- 
sive (M < 3.3 x 10 3 M Q ) than CI. 
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